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We derive an analytic expression of the effective potential at finite temperature (T) and 
chemical potential (/x) in the strong-coupling lattice QCD for color SU(3) including next- 
to-next-to-leading order (NNLO) effects in the strong coupling expansion. NNLO effective 
action terms are systematically evaluated in the leading order of the large dimensional (1/d) 
expansion, and are found to come from some types of connected two plaquette configurations. 
We apply the extended Hubbard- Stratonovich transformation and a gluonic dressed fermion 
technique to the effective action, and obtain the effective potential as a function of T, \i, and 
two order parameters; chiral condensate and a vector potential field. The next-to-leading 
order (NLO) and NNLO effects result in modifications of the wave function renormalization 
factor, quark mass and chemical potential. We find that T CjM =o and /i C) T=o are similar to 
the NLO results, whereas the position of the critical point is sensitive to NNLO corrections. 

§1. Introduction 

Understanding the Quantum Chromodynamics (QCD) phase diagram is one of 
the most interesting problems in quark and hadron physics. In the present Rclativis- 
tic Heavy-Ion Collider (RHIC) experiments, phase transition to strongly coupled 
matter consisting of quarks and gluons seems to be observedP The QCD phase 
transition observed at RHIC takes place at almost zero baryon density, where the 
predictions of lattice QCD Monte-Carlo (MC) simulations are reliable. The phase 
transition of compressed baryonic matter will be probed in the future experiments 
in FAIR, J-PARC and low energy programs at RHIC. At finite baryon densities, the 
lattice MC simulations are difficult due to the complex fermion determinant P In 
order to discuss the whole shape of the phase diagram, it is necessary to invoke some 
approximations in QCD or to apply effective models. The strong-coupling lattice 
QCD (SC-LQCD) is one of the most instructive approximations to investigate the 
phase structure at finite temperature T and chemical potential \i. 

SC-LQCD was first applied to the pure Yang-Mills theory. Wilson showed that 
the Wilson loop would follow the area law at strong coupling,® and the Creutz 
demonstrated that the lattice MC simulation can connect strong coupling and weak 
coupling (perturbative) expressions of the string tension.® The behavior of the 
string tension as a function of the inverse coupling, j3 = 2N c /g 2 , is well under- 
stood in the strong coupling and character expansions, which were developed by 
Minister.® Chiral symmetry in SC-LQCD also has been long studied from 1980s. 
Basic formulations has been developed based on the staggered,® >® Wilson®® and 
naive® fermions. The domain-walP® and the overlap^ 1 fermion provide modern 
formulation of the lattice chiral symmetry, and some SC-LQCD based investigations 
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are found in Refs. fl~2j) JT3]) (domain- wall) and fT4" |) . [T5]) (overlap). In the strong cou- 
pling limit (SCL), chiral symmetry is spontaneously broken in vacuumj^^'^^QS 
and restored at high T and/or large ^QSi-EE} 

Since the SC-LQCD is based on the same formulation as the lattice MC sim- 
ulations, its results should be consistent with the MC results such as hadron mass 
spectrumP^~® as long as the applied approximations are valid. The phase dia- 
gram structure has been predicted in the strong coupling limit ,E7J) anc j ^ j g 
recently confirmed qualitatively in MC simulations® based on the monomer-dimer- 
polymer (MDP) formalism.® In order to make a step forward towards the true 
phase diagram, it is necessary to develop the formalism to include the plaquette 
effects both in lattice MC simulations at finite [i and SC-LQCD. There exists some 
SC-LQCD works including NLO effects, i.e. one plaquette contributions, on the 
hadron masses® ~® and the phase diagram EDiEDJSSJ i n our previous work on the 
NLO SC-LQCD,® we find that the phase dia gram evolves to an empirical shape 
with increasing ft = 6/g 2 , while the critical temperature at zero chemical potential 
T c ,[i=o is larger than the MC results. The later observation suggests that we need 
to evaluate the next-to-next-to-leading order (NNLO) effects on the phase diagram, 
which have never been investigated before. 

In this paper, we derive an analytic expression of the effective potential including 
NNLO effects at finite T and fi, and investigate NNLO contributions to the phase 
diagram. We adopt one species of unrooted staggered fermion corresponding to 
Nf = 4 in the continuum region. Since the flavor dependence of the phase boundary 
has been shown to be moderate,®'®'® the present results could be valuable for 
the understanding of the phase diagram with Nf = 2 + 1. Effective action terms 
from one and two plaquette configurations are obtained by integrating out spatial 
link variables. We apply the extended Hubbard-Stratonovich transformation®'® 
to bosonize fermion interaction terms. With the l/<? 4 corrections, we encounter 
those terms containing the next-to-nearest neighbor (NNN) interaction, which can 
be evaluated by introducing a gluonic dressed fermion. The effective potential is 
obtained as a function of T, /x and two order parameters: the chiral condensate 
and the quark number density. We determine the equilibrium from the stationary 
condition of the effective potential with respect to the auxiliary fields, and study 
the properties of QCD phase diagram. MC studies based on one species of unrooted 
staggered fermions have been carried out extensively around ft ~ sPD-BSD For the 
comparison between SC-LQCD results and those in the MC simulations, we discuss 
the results in the region ft < 6. As shown later, the comparison of NNLO and 
NLO results suggests that ft values under consideration are in conversion radius. We 
compare the critical temperature at [i = (T Cj u=o) with the NLO and MC results. 
We also study the evolution of the phase diagram and the critical point, as well as 
the possibility to have partially chiral restored matter, which has been suggested in 
NLO SC-LQCD.®'® 

The QCD phase transition has another aspect of the deconfinement transi- 
tion. In order to discuss the deconfinement transition, it is necessary to include 
the Polyakov loop effects as discussed in the framework of SC-LQCD®'®'® and 
in the Nambu-Jona-Lasinio model with the Polyakov loop (PNJL) E3Ji,[54j j n 
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study, we concentrate on the chiral phase transition, and the simultaneous descrip- 
tion of chiral and deconfinement transitions will be reported elsewhere P3 

On the phase transition study with unrooted staggered fermion, we still have 
some debate. In the strong coupling limit, the effective action with one species of 
unrooted staggered fermion is in the same universality class as the three-dimensional 
0(2) spin models, and the phase transition is the second order In the continuum 
region, one species of unrooted staggered fermion corresponds to four flavors (or 
tastes). The phase transition at // = is expected to be the first order for Nf > 3 
from the anomaly argument.^ The axial anomaly may cancel from the doublers on 
the lattice,^ but it is discussed that the conventional staggered anomaly appears in 
the taste-singlet PCAC relationPS 1 In numerical MC simulations pHJ the transition 
at pi = is shown to be the first order in the continuum region. While we have these 
complications, the staggered fermions has a merit that it is simple. As a result, it is 
fast in numerical simulations, and it is possible to carry out the analytic calculations 
in the mean field approximations when we adopt the unrooted staggered fermion. 

This paper is organized as follows. In <j2l we derive the effective potential with 
1/g corrections. In we show the calculated results of the effective potential and 
the phase structure. We summarize our work in $3} All through this paper, we 
use the lattice unit a = 1, and physical values are shown in dimensionless values 
normalized by the lattice spacing a. 



§2. Effective Potential in NNLO SC-LQCD 



In this section, we derive the effective potential in SC-LQCD with NNLO cor- 
rections (1/g 4 ) at finite temperature T and chemical potential \i. We start from 
the lattice QCD action with one species of unrooted staggered fermion (rif = 1) 
for color SU(N C ) (§ 12. which corresponds to Nf = Arif = 4 in the continuum 
limit. After a short review on SCL and NLO effective action (§ 12. 2p . we derive the 
NNLO effective action (§ 12. 3p . Interaction terms are reduced into a bilinear form of 
fermions through the extended Hubbard-Stratonovich transformation (§ I2.4p and by 
introducing a gluonic dressed fermion (§ 12. 5p . We obtain an analytic expression of 
the effective potential in 



2.1. Lattice QCD action 

The lattice QCD action and the partition function with one species of staggered 
fermion for color SU (N c ) are given as follows 



^lqcd = / T>\x,xM e" SLQCD = J V[ X ,x,Uu] e"^'-^-^ , 



~ 2 y~iy~i%> xxUj :X x x+ ~j x x+ jU] tX xx 

X j = l 



(2-1) 
(2-2) 

(2-3) 
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^ = 4EE[ tr u ^ + tr <J - 3 E E [ tr ^-m + tr u k J > ( 2 - 4 ) 

a j>0 y a- fc>j>0 

where x(x)j m 0i Uo tX (Uj jX ), and ?7oj,a: (Ujk,x) denote the quark (anti-quark) field, the 
bare quark mass, temporal (spatial) link variable, and temporal (spatial) plaquette, 
respectively. The spinor structure is compressed into the staggered phase factor 
Vj,x = (— l) 2-0 " 1 ' rX i~ 1 _EJ.ED.I12J We have introduced two types of mesonic composites 

V+ =Xxe^U 0txXx+ o , V- = X x+ de-^Ul xX x , M x = XxXx , (2-5) 

which appear in the effective action discussed later. Quark chemical potential fi on 
the lattice is introduced as a weight of the temporal hopping in p^JSJ By using 
a 75-related factor e x = (— l) xo_l a staggered chiral transformation is given as 

Xx e ietx Xx®^^^^ The lattice kinetic action Sp is invariant under this chiral 
transformation in the chiral limit tuq — > 0. Throughout this paper, we consider the 
case of color SU (N c = 3) in 3+1 dimension {d = 3) spacetime. Temporal and spatial 
lattice sizes are denoted as N T and L, respectively. While T = 1/N T takes discrete 
values, we consider T as a continuous valued temperature. We take account of finite T 
effects by imposing periodic and anti-periodic boundary conditions on link variables 
and quark fields, respectively. We take the static and diagonalized gauge (called 
Polyakov gauge) for temporal link variables with respect for the periodicity.^ In 
these setups, we evaluate the effective potential Teg = — log [^lqcd] / (N T L d ) based 
on the strong coupling expansion. 

2.2. Strong coupling expansion and SCL and NLO effective action 

In a finite T treatment of SC-LQCD, we first derive an effective action by inte- 
grating out the spatial links Uj, and the temporal link variable Uq is evaluated later. 
The Uj integral can be exactly performed in each order of the 1/g 2 expansion by 
utilizing the link integral formulae 

/ dU U ab U\ d = ±-5 ad 5 bc , / dU U ab U cd ... U ef = jLe ac ... e £b ,.. f , (2-6) 

and so on. Then an isotropic hopping structure of hadronic composites emerges, 
and enables us to simplify the effective action by utilizing another expansion: 1/d 
expansion, which is explained in the next subsection. 

In this work, we consider the effective action including the SCL (l/g°), NLO 
(1/g 2 ) and NNLO (1/g 4 ) terms in the strong coupling expansion, while we keep only 
the leading order terms (l/d°) in the 1/d expansion. The effective action is defined 
as 

e -S eS (x,X,Uo) = j vu . e -5 LQC D = e -4 T) J VUj e - S F ) - s G 



=e sP-4 s L ( e -s G ) . (2 . 7) 
We have defined an expectation value as 

(O) J VU, 0[Uj] e~^ , Z« = / WJ e- Si;) = ^ ■ M 



^SCL 
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Fig. 1. Diagrams contributing to the SCL and NLO effective actions. Quarks (anti-quarks) are 
shown in open (filled) circles, and spatial (temporal) link variables are represented by solid 
(dashed) lines. Note that we should also take account of hermite conjugate contributions for 
NLO diagrams. 



The factor l/Zg^ L ensures the normalization property (1) = 1. In order to sys- 
tematically evaluate the effective action terms in each order of 1/g 2 , the cumulant 
(or coupled cluster) expansion is indispensable. It is well known that the expecta- 
tion value of the exponential form operator with a small factor (i.e. 1/g 2 ) can be 
evaluated by using the cumulant expansion® 1 



So 



n=0 



111 



[s G n ) 



exp 



n=l 



(-1)' 



III 



[Sg) c 



The correlation part in the connected diagram contributions is shown by the bracket 
(• • • ) c , and is called a cumulant, e.g. (S G ) C = (Sq) — {Sg} 2 - By substituting Eq. (|2-9p 
into Eq. (j2-7|) . we find that the effective action is obtained as 



S c s =5scl - 



n=l 



111 



-Sscl + 4S NLO + ^Snnlo + 0(1/ g 6 , 1/Vd) . 



(2-10) 



The n-th term in the sum is proportional to l/g 2n , and we can identify n = 1 and 
n = 2 terms as NLO and NNLO effective action, Z\5nlo and ^SnnlO) and Sscl 
shows the SCL effective action. 

The SCL and NLO effective actions are known asp^'® 



Sscl = 
^Snlo 



4 r) + 4 S) 



-s 



(r) 



1 



(2-11) 



x,j>0 



G)c 



4N^g 2 



-2 x 



£K 4 



V- , + V + ,TC 



x,j>0 
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X+J 



X+J 



16iV c V 



(2-12) 



x,k>j>0 



The factor and sign of each term is summarized in Table [H The factor "2" in 
the last line accounts for the hermite conjugate contribution. In Fig. (TJ we show 
the diagrams contributing to the SCL and NLO effective actions. The first three 
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diagrams represent the temporal hopping (Sp ) terms and the spatial meson hopping 
(MM) terms in iSscl- The fourth and fifth diagrams show VV and MMMM terms 
in Z\Snlo 

2.3. NNLO effective action 

From Eq. ()2- lOj) . we find that the NNLO effective action is given as the cumulant 



where Up = txU^ v>x denotes the trace in the color space of a plaquette, or its conju- 
gate. In the case where the two plaquettes (P and P') do not have any common spa- 
tial link variables, the average of the product is factorized as (UpUpi) = (Up)(Up>), 
which cancels with the second term. Thus only the connected diagrams can con- 
tribute to the NNLO effective action. In Fig. [21 we show the two plaquette configu- 
rations, where the plaquettes share at least one spatial link. 

For a given two plaquette configuration, we consider the leading order terms 
in the 1/d expansion.® The sum over spatial directions in Eq. (|2-lip would 
give rise to a factor d due to the spatial isotropy. Provided that the meson hopping 
term V . M X M^.,~- in the SCL effective action stays finite at large d, the quark field 

(x,x) should scale as <i -1 / 4 . Then a factor d from the sum over spatial directions 
and a factor 1/d from the four quark fields cancel, and the mesonic hopping term 
Y2j M x M x+ j is found to be O(l/d ), which is the leading order in the 1/d expansion. 
This also applies to the NNLO diagrams. For example, we consider the product of 
the same temporal plaquette, NNLO-TT1 in Fig. [2j We will have a factor d for the 
sum over the temporal plaquette Uoj tX , and four quark fields give rise to a factor 
1/d. Thus the diagram NNLO-TT1 in Fig. [2] is in the leading order 0(l/d°) in the 
1/d expansion. Diagrams with more quarks for the same plaquette configuration are 
suppressed as 0(l/y/d) for N c > 3. This is called the systematic large dimensional 
or 1/d expansion. 

Since we concentrate on the leading order terms in the 1/d expansion, we find 
that it is sufficient to consider the diagrams shown in Fig. [21 each of which includes 
the minimal number of quark fields for each plaquette configuration. We omit two 
plaquette configurations with Up = Up,, which lead to a constant in the effective 
action in the leading order of the 1/d expansion, 0(1/ (P). By utilizing the group 
integral formulae Eq. (|2-6p . these contributions lead to the following effective action 
terms 



of 52 



^Snnlo = - \{sh) c = -\ [(sl) - (S G ) 2 } 




P,P' 



(2-13) 




(2-14) 
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Fig. 2. Diagrams contributing to the NNLO effective action. Definitions of the symbols and 
lines are the same as those in Fig. [1] The top, middle and bottom rows represent the NNLO 
diagrams composed of the temporal-temporal, spatial-spatial, and temporal-spatial plaquette 
configurations, respectively. 
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UN 7 Q A E MxM^M^M^M^M^ , (2-15) 

c!J fc>0, \j\jtk, \l\^\j\,k 



1 



16iV c V 



x,j>0,\k\^j 
^ ^x+j+k^x+k 



v+v- , + v + ,K" + v + ..y- ft + y+ .y- . , 

a;— 0+j x— x— x— 0+j 



(2-16) 



The factor "2" in the first line of Z\5^ NLO accounts for the hermite conjugate con- 
tribution. We find that two new mesonic composites appear 



w+ = w 2 "rn.,r t ,,.. oV ,. 



o , W- = X x+2 ^Ul x+b Ul xXx . (2-17) 



These composite connect the quark fields in the next-to-nearest neighboring (NNN) 
temporal sites. In Fig. [5J the top, second and bottom rows represent AS^^ LO , 
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Table I. The factors and sign which come from, (A) a global NNLO factor l/g°, 1/g 2 or l/g 4 , (B) 
a symmetric factor (Sym.) 1/2 for the product of two same plaquettes, (C) product of factors 
— 1/2 and 1/2 for the forward (XUx) an d backward (XU x) fermionic elements (F.E.), (D) a 
product of the staggered factors T]j tX , and (E) one- link integrals. Each entry in the left column 
represents a diagram shown in Fig. [T] and [2] 





(A) Global (B) Sym. (C) F.E. (D) (E) U 3 


Total 


SCL-S 


l/g° 1 -(1/2) 2 +1 1/JV C 


-1/4N 


NLO-T 
NLO-S 


1/g 2 1 -(1/2) 2 -1 1/N 2 C 
l/g 2 1 + (l/2) 4 -1 1/Nt 


+ l/AN 2 c g 2 
-l/lQN 4 c g 2 


NNLO-TT1 

NNLO-TT2 

NNLO-SS1 

NNLO-SS2 

NNLO-TS1.2 


1/ff 4 1/2 -(1/2) 2 -1 1/iV 2 
l/g 4 1 -(1/2) 2 +1 l/N? 
l/g 4 1/2 +(l/2) 4 -1 l/Nt 
l/g 4 1 -(1/2) 6 +1 l/Nl 
l/g 4 1 +(l/2) 4 +1 l/N? 


+ l/8N 2 c g 4 
-l/4N?g 4 
-l/32N 4 g 4 
-1/64JV C V 
+l/16iV c V 



^^nnlo an< ^ ^^nnlo' respectively. Note that for diagrams which include quarks, 
we also consider the hermite conjugate of those diagrams. 

We summarize the factor and sign in SCL, NLO and NNLO in Table[H The factor 
and sign of each term is a product of (A) a global NNLO factor l/g°, 1/g 2 or 1/g , 
(B) a symmetric factor 1/2 for the product of two same plaquettes, (C) product 
of factors —1/2 and 1/2 for the forward {xUx) an d backward {j(U x) fermionic 
elements, (D) product of the staggered factors m x , and (E) a factor from the one- 
link integral. Expressing the effective action as a function of color singlet composites 
M, V and W, we generally obtain a factor "—1" resulting from an odd number 
times of grassmann number exchanges. This sign factor is taken care in the global 
factor and shown in the column (A) of Table HI For a product of different plaquettes 
(P, P'), there are two combinations (P, P') and (P' , P) in the sum of Eq. (|2-13|) . and 
they give rise to a factor 2. With the factor — l/2g 4 found in Eq. (|2-13p . we get a 
global factor 1/g 4 . For a product of the same plaquettes P = P' in NNLO-TT1 and 
NNLO-SS1, we need a symmetric factor 1/2. The factor from the one-link integral 
is 1/N*, where k is the number of integrated out spatial links. Also in the case to 
use the second formula in Eq. (|2-6p . we have 1/N C in the present configurations. For 
example, in NNLO-TT1, we have a factor l/(iV c !) 2 from the spatial link integrals, 
and the following combination of temporal link variables remains on the link (0, x+j) 

e aibl ... Cl e a2b2 ... C2 U b ' h2 , ■ ■ ■ XJ C4C2 , = (N c - 1)! x Ul a2a i . (2-18) 

We also find another (N c — 1)! appears on the link (0,x), and the factor from the 
one-link integral is found to be 1/N 2 . In a similar way, 1/N 4 appears for NNLO-SS1. 
Note that NLO-S and NNLO-SS1 action terms are hermite, hence one find a factor 
2 in MM MM terms of Eq. (l2T2"1) and (I2d5]) . 

We explain the factor and sign for NNLO-TS1, as an example. There are two 
forward ((k, x) and (j, x + k)) and two backward ((fe, x + j) and (j, x + 0)) fermionic 
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elements, and we have a factor +(l/2) 4 . Staggered factors for these elements are 
Vk,x, V jjX+ k = dr] jtX , Vk, x +'j = ~ a Vk,x and 7? J>+Q = -rj j>x , where a = 1 (-1) for k > j 
{k < j). The product of these staggered factors results in +1. Five spatial links are 
integrated out, and give a factor Together with the global factor 1/g 4 , we find 

the coefficient +l/16iV 6 5 9 4 for NNL0-TS1. 

We combine those terms in NLO and NNLO having the same composites. Now 
we have obtained the effective action S e g in Eq. (|2-l(jp up to 0(l/g 6 , 

S e & = Sscl + AS T + AS S + AS TT + AS SS + AS TS (2- 19a) 

= \ - KT) - \d ? [MM]j ' x (2 ' 19b) 

x x,j>0 

+ \Td^ [v+v ~ + v ~ v+kx (2 ' 19c) 

x,j>0 

~ldd~T) S [MMMM] jk , x (2-19d) 

^ ' x,j>0,k>0,k=ij 

\ w+w ~ + W-W + ]j, x (2-19e) 

x,j>0 

v /v y i,j>0, |fc|>0, |Z|>0 

x,j>0,\k\^j 

We have introduced a short-hand notation 

L4i% x = A x B x+] , [ABCD] jkjX = A x B x+ f x+ , +k D x+k . (2-20) 

We obtain Eq. Q2-19gP by shifting x to x + in the third and fourth V + V~ terms 
in Eq. ([2TS]) . NLO and NNLO effective action terms, /IS^A" = r, s, rr, ss, ts), 
correspond to Eqs. (|2-19c|) -( |2-19gp . The coefficients are defined as 

8 8 - W-l)(d-2) , - d{d ~ l) (2-22) 

2.4. Extended Hubbard- Stratonovich transformation 

The effective action S^ff derived in the previous subsection contains several types 
of composite products, which include many fermion fields. In order to perform the 
grassmann integral, it is much more convenient to reduce these higher order terms to 
the spatially local and bilinear form in the fermion fields through the bosonization 
procedure, so-called the Hubbard-Stratonovich (HS) transformation. The standard 
HS transformation is applicable to the product of composites of similar kind, while 
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Table II. The application of the Extended Hubbard- Stratonovich transformation to effective action 
terms. 





Q 






AS TT 


(|2-19e|) 


/W2d 


w+ 




AS SS 


(|2-19f|) 


p ss /4d(d - l)(d - 2) 




[MMM%,, 


AS TS 


<2-19g) 


P TS /8d(d - 1) 


-[v+v-]^ 




AS T 
AS 8 


H2-29I) 
(12-301) 


(/3 T +/3 TS ^T S )/4d 
(ft + 2p as i> sa )/2d{d - 1) 


-v+ 

[MM],-,, 





NLO and NNLO terms contain the product of different types, such as V + V~ . In 
order to treat these terms, we apply the extended Hubbard-Stratonovich (EHS) 
transformation,^'^ where two auxiliary fields (<£>, tfi) are introduced simultaneously 



d v d4>e- a {^- {A+B)/2]2+ ^- i{A - B)l2]2 } +aAB 

e -a{ v 2 ~(A+B) v +^ 2 -i(A-B)<f>} 
^—a^-ipip—Aip—ipB} 

stationary 



stationary 



(2-23) 
(2-24) 



where a is an arbitrary positive real constant. In the second line, we approximate 
the integral by the integrand at the stationary values, ip = {A + B)/2 and <j> = 
i(A — B}/2 (the saddle point approximation). This procedure is denoted by 
and we use this notation in later discussions. In the last line Eq. ()2-24p . we have 
performed the transformation tp = ip + id) and tp = ip — i<f>. When A = B is satisfied, 
Eq. (|2-23p is equivalent to the standard HS transformation. In the case where both 
(A) and (B) are real, the stationary value of <f> is pure imaginary. Thus we replace 
(f> — > iu in Eq. ()2-23p . then the stationary value of uj is real 



nAB 



-a{ip 2 -(A+B)ip-u 2 +(A-B)uj} 



stationary 



(2-25) 



We first apply the EHS transformation in the complex representation, Eq. (|2-24p . 
to NNLO terms in Snnlo, Eqs. (|27T3el) . $2lM and ( |2-19g| ). By substituting (a, A, B) 
in Eq. (|2-24p as shown in Table HH we obtain 



AS S 



x,j>0 

~N T L d /3 TT i; TT ilj TT 

~4d(d-l)(d-2) 

~N T L d p ss 4> ss il: ss - 



(2-26) 



x,j>0, \k\>0, | Z | > 
\k\?j,\l\&,\l¥\k\ 

- ^ ss [MM], 

x,j>0 



i> ss ^ ss - [MM] jjX+i ^ ss - $ ss [MMMM] jKl 

fiss^ss 



d(d - 1) 



[MMMM] jkyl 



x,j>0,k>0,k^j 



(2-27) 
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^ "MfTTTT E firsts + [T^k^ 
1 ' x,j>0,\k\^j 

-i s ([M] jii+£ + [MM]. 

(2-28) 

We have introduced the auxiliary fields (^k^k) in (If = tt,ts,ss). In 

Eq. (|2-26p - (|2-28p . the symbol "~" represents that constant and isotropic values 
are assumed for auxiliary fields (iPk^k), and we use this notation in later dis- 
cussions. Note that the EHS transformation is applied so that the auxiliary fields 
(tpss^ssjiprsj^rs) have no chiral charge in Eqs. (|2-27p and (|2-28p . For example, 
in Eq. (|2-27p . we have decomposed the six-meson term {MM MM MM) into MM 
and MMMM mesonic terms, and introduced auxiliary fields for their composites, 
(ipss^ss) = ((MM), (MMMM)). 

NNLO contributions in Eqs. (p 7 ^ and (flZTZHj) generate terms, [V + V~ + 
and [MMMM]. They can be absorbed into the NLO contributions, Eqs. (|2-19cl) 
and (pTOp . We put [V + V~ + V~V+] and [MMMM] terms in NLO and NNLO 
together and bosonize as 

A ~ T = P T + firsts ^ + y-y+^ 

x,j>0 

0' r - - i 

x,j>0 

^N T L d ^-^ T + ^ ^ T V+ - AV-) , (2-29) 

x 

A ~ s = _ ^ [MMMM] jk , x 



[Ms - [MM] j>x ip s - $ s [MM] j 



^ ' x,j>0,k>0,k^j 

~N T L d ^ s t(> s -^Y(^s + i)[M% , (2-30) 

x,j>0 

where 0' T = T + TS ipr S and 0' s = S + 20 ss ijj ss . We find that V + V~ and MMMM 
terms shifts the coefficients of SCL effective action terms via EHS transformations. 
In previous NLO investigations JSJJHUJ we have developed the method to express 
V + V~ and MMMM effects as modifications of the wave function renormalization 
factor, quark mass and chemical potential. In the present investigation, we need 
to manipulate the next-to-nearest neighboring interaction come from effects in 
Eq. (|2-26p before utilizing those techniques developed in NLO. The effects are 
evaluated in the next subsection. 
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We introduce the chiral condensate auxiliary field a by using the standard HS 
transformation. We combine the mesonic hopping [MM] terms in 5 SC l (Eq. (|2-19bjl ). 
NNLO contributions (Eqs. §FT7\ and and NLO contributions with NNLO 

effects (Eq. (|2-30p ). Then we obtain the modified mesonic hopping terms, which can 
be bosonized as 



f£M x V xy M y ttb' a Y^ 



1 



cr x V xy a y + a x V xv Mi 



x ' xy lv± y 



b' 



^N T L d ^o 2 + b' a oY,M x , 

X 

b' a =b a + 2 [P s s4> ss + P TS 4> TS + (p s + 2 / 3 ss ^ ss )(^ s + $„)] 

V *y =^1 E^+5',2/ + 6 x-j,y) ■ 



(2-31) 

(2-32) 
(2-33) 



We have obtained the additional mass term b' a a £V M x . Combined with the current 
quark mass mo, we obtain a constituent quark mass as 



b' a a + m 



(2-34) 



Now we have the NNLO effective action in the spatially local and bilinear form 
in the fermion fields. We collect all the remaining terms from Sscl (Eq. (|2-19bD ). 
NNLO terms (Eqs. Q£2HD , (|2^7|l and (J23SD), NLO terms (Eqs. (IMj) and (fCTD ). 
and the chiral condensate term (Eq. (|2-3ip ) 



c _c(f) i cW 



C =1 + (/3 T + /3 ts Vts)?/V , C = 1 + (J3 T + P TS lpTs)i>T 



(2-35) 
(2-36) 

(2-37) 



S. 



(X) 
cff 



--N T L a 



PttIPtt^tt + Pss^Pss^ss + ^rsVVsVVs + ^(Pt + Pts^ts^t^t 



+ l((3 s + 2p ss ip ss )ij; s ip s + \b' a a 2 



(2-38) 



where and S*^ represent those effective action terms with and without fermions, 
respectively. The SCL and NLO effective action contains V x , V~ and M x , and their 
corresponding effective potentials are known.® In order to obtain the effective 
potential from the NNLO effective action containing W x and W~ in addition to 
Vj~, V~ and M x , we discuss how to treat W+, W~ terms in the next subsection. 

2.5. Gluonic dressed fermion 

The fermionic NNLO effective action S)^ in Eq. (flZ^o]) is in the spatially local 
and bilinear form of x anci X> then it is possible to obtain the fermion determinant 
analytically in the form of Matsubara product. In order to evaluate the Matsubara 
product, we need to obtain the Matsubara frequency which gives zero determinant. 
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Because of the coupling of the next-to-nearest neighboring (NNN) temporal sites via 
W + and W~, the fermion matrix becomes pentadiagonal rather than tridiagonal 
and it is not easy to obtain the solution. In the present case, NNN terms are in the 
0(l/g 4 ), then we can absorb them by introducing a gluonic dressed fermion 

X ' x = Xx + Ae»U ,x X x+ 6 > Xx = Xx + A^e^U^ > ( 2 ' 39 ) 

where A = 0(l/g 4 ). Mesonic composites are represented in this dressed fermion as 

M x =M' X - AV'+ - AY- + 0(l/g 8 ) , (2-40) 

V+ =V'+ - AM' x+6 - AW' X + + 0{l/g 8 ) , (2-41) 

V- =V<- - AM' x+d - AW- + 0{l/g 8 ) , (2-42) 

W+ =W' X + + 0(l/g 8 ) , W~ = W'~ + 0(l/g 8 ) , (2-43) 

where primed composites denote those with dressed fermions, such as M' x = x' x x' x - 
In terms of this dressed fermion, the fermionic term S^g is rewritten as 



X 



\{CA-CA) 



ML 



~\ E [( 2 ^rV>rr + C A)W X + + (2/3 TT ^ rr - CA)W X ~] + 0(l/g 8 ) . (2-44) 

x 

We find that NNN terms in the primed fields cancel with each other when we choose 
A = —2(3 TT ip TT /C and A = 2f3 TT ijj TT /C. Then the NNLO fermionic effective action 
is represented in V x + , V x ~ and M' x terms. In the later discussion, we regard x' an d 
X 1 as quarks and anti-quarks, and omit the prime for the mesonic composites. The 
fermionic effective action now reads 

S iP =\ E i Z - V * + ~ Z + V *] + E m 'i M * + °(V5 6 ) , (2-45) 

X X 

Z_ =C - 2m q A = C + 4m' q (3 TT Tp TT + 0(1//) , 

Z + =C + 2m q A = C + 4m' q P TT ^ TT + 0(1/ g 6 ) , (2-46) 
1 

2 1 

Comparing Eq. ()2-45p with (|2-36|) . we find that the effects of are expressed as 
modifications of the constituent quark mass and the coefficient of through the 
gluonic dressed fermion. It should be noted that the Jacobian between x) and 
(x')X') deviates from unity by 0(1/ g iNrNc ), and we can ignore its effects in NNLO 
SC-LQCD. As a result, NNLO effective action is represented in Vj~,V~ and M x 
terms as in the case of bosonized NLO effective action.®'^ 

2.6. Effective potential 

In the discussion by the previous subsection, we find that the NLO and NNLO 
corrections lead to the coefficient modification of V + , V~ and M in the fermionic 



m' q =m q --(CA- OA) = m q - (3 TT (^ TT + Vw) + 0(l/g 6 ) . (2-47) 
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effective action. The meaning of this modification would be more clearly understood 
in the following representation of the fermionic effective action 



2 



s%> =\ E [ z - v * - z +v*] + E < M > 

; x 

E \ - + E 

- x iV x 

=Z x E * x G xy(-™-<i> M> T )Xy , 

xy 



rh q M x 



Wlq&xy i 



111 



Z x = x/Z^ZZ , m q = —1 , p, = n - S/J, = \i - log y/Z+jZ- 



(2-48) 

(2-49) 
(2-50) 



In this form of the effective action, we easily find that NNLO effects result in the 
modification of the wave function renormalization factor Z x , quark mass rh q and 
chemical potential fl by using the technique developed in NLO investigations jMJJSS) 
We carry out the Grassmann integral over the quark fields and the temporal link 
integral in a standard wayP^^^ First, we perform the Fourier transformation in 
the temporal coordinate, and obtain the product in the frequency by the Grassmann 
integral over the quark fields. Second, we evaluate the product in the frequency by 
using the Matsubara method. Finally, we carry out the temporal link integral by 
using the Haar measure. The fermionic part in the effective potential jF„ is given as 



N T L d 



log 



J V[ X ,x,U ]e~ s ^ 



V q (m q ;fl,T)-N c logZ x , (2-51) 



V q (fh q ; fl,T) = -Tlog 



smh[(N c + l)E q (m q )/T] 
smh[E q (rh q )/T] 



+ 2cosh(N c fl/T) 



(2-52) 



where the number of temporal sites is replaced as N T = 1/T, and E q {m q ) = 
arcsinh (rh q ) is regarded as the quark excitation energy. 

The latter point is understood from the expression of the temporal quark hopping 
matrix determinant 



-V q L d /T 



Jv[U ]det XtC G^(m q ;fi,T) (2-53) 
[] / ^odet c \e E <>/ T (l +U e-( E «-M T ) (l +aJ e -( J5 ,+/i)/T> 



(2-54) 

where ^o( x ) = Y\ T £^o( x j f). We can evaluate this temporal link integral exactly, for 
example, in the Polyakov gauged 



(2-55) 
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We find the Boltzmann factor exp [—(E q ± ji)/T] appears, and the integrand is noth- 
ing but the thermal quark contribution to the partition function if Uq is close to unity. 
Note that the determinant is for the spacetime and color in Eq. (|2-53|) . and for the 
color in Eq. (j2-54p . 

We obtain the NNLO effective potential as follows. 



fi,T)+J*x\$) , (2-56) 
J-W =SSP/N T L d 

+ Ptt^tAtt + Pssi>ss4>ss + ^Pts^ts^ts , (2"57) 

m! t 

rh q =— 1 , m = b a a + m - (3 tt (i/j tt + ip TT ) , (2-58) 
z x 

Z + =1 + {p T + P TS ip TS )ij; T + 4/3 TT m q $ TT , (2-59) 

Z_ =1 + (/3 T + p TS tjj TS )^ T + 4/3 rr m^rr , (2-60) 

b' a =b a + (p, + 2P 8a $ M )(rl> a + $ a ) + 2p ss ^ ss + 2/W rfl , (2-61) 

where <P = (a, {i/jk^k', K = t,s, tt, ts, ss}). Coupling constants (b a ,f3 T , f3 s , (3 TT , j3 TS , j3 s 
are defined in Eqs. (EF2TT) and (l2 T 22]l . 

2.7. Stationary conditions 

The auxiliary fields introduced during the bosonization procedure have to satisfy 
the stationary condition, dJ- c ^jd^> = 0. Within the constant field approximation, the 
stationary values of auxiliary fields are related with each other, then we replace the 
auxiliary fields other than (u, ip T and ip T ) with their equilibrium values as functions 
of (cr, Vv and ip T ). 

We can solve the stationary conditions for all auxiliary fields exactly. In solving 
the stationary condition, we first note that the quark free energy is a function of m' q , 
Z + and Z_, as we can guess from Eq. (|2-45p . Then the variation J- q is obtained as 

dF q _ dF q dm' q dF q 8Z+ dF q dZ- 

d<P dm' q d<P dZ + d<P dZ_ d<P { ' 

Since Z + and Z_ contain m' , the above derivative dJ- q /dm' q also contains the deriva- 
tive via Z±, (dJ 7 q /dZ±)(dZ±/dm' q ). Substituting a,ip s , , ip s ,'ip T ,ip T , in Eq. (|2-62[) . the 
stationary conditions for these fields read, 

By using these relations, the stationary conditions for NNLO auxiliary fields (V> TT , 
ip TT , ip TS , ip TS , ip ss , ip ss ) are obtained as 

■tjj TT = 2m' ijj T — o , ijj TT = 2m'ijj T — a , (2-64) 
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4>ts = VvVv , 4>ts = 2a 2 , $ ss = a 2 , tp ss = a 4 . (2-65) 
From Eq. (|2-64p . we obtain function of (a, ip T ,ip T ) 



j _ h' a a + m + 2(3 TT a 
l + 2(3 TT {A + A) 



m q = - — — — — — — ■ — . (z-DDj 



Now all the stationary values of auxiliary fields other than (<r, f/Vi Vv) are obtained 
as functions of these three fields. 

We shall now proceed to erase one more auxiliary field. Three derivatives in 
Eq. (|2-62p . dT q jdm! q and 8J r q /dZ±, are not independent, but are related to the 
two derivatives, dV q /dm q and dV q /dfi. The stationary conditions for tp T and ip T in 
Eq. (|2-63p are rewritten in terms of these derivatives as 

Z-A=^j^~^ + N c (2-67) 
Z x om q ofi 

m' dVa dVn 

We solve the coupled equations, Eqs. ()2-67p and (|2-68p . for dV q /drh q and dV q /djl 



dV, 



q - - u T (l - Am' q (3 TT a) , (2-69) 



=<p T + /^VvVv + 2p TT m' q (i> TT i) T + ArA) ~ N c . (2-70) 

We have parameterized the auxiliary fields Vv ; Vv as tjj T = (p T + uj t and tp T = tp T — 
u T . This corresponds to applying the EHS transformation in Eq. (|2-25p . where we 
substitute (f> — > iuj. Equation (|2-69j) implies that oj t is related to the quark number 
density, p q = —dF c s/dp = —dV q /djl. It is also possible to solve the stationary 
condition for a in Eq. (|2-63p for dV q /dm q . 

~^T^h = ~ <T + WtMttA + VVrVv) , (2-71) 

Two equations Eqs. (|2-70|) and (|2-7ip have to be consistent, then we find the relation 
between a, lj t and ip T 

( p T + P> T ( (p 2-u;2) + m q a-N c = . (2-72) 

By substituting Eq. (|2-66p . we obtain a cubic equation for <p T 

A(3' T Ard + (P'r + 4/3 TT + 2/3 rs <r 2 )<^ + (1 - 4iV c /3 rr - 4/3;/3 rr ^) v?T 

-N c + a[(b a + 2(3 TT + 2(3 s a 2 + 6^ ss a 4 - 2^ TS u? T )a + m ] - fcu 2 = , (2-73) 

and its solution ip T = Lp T (a,uj T ) is obtained analytically. 
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We have obtained the final expression of the effective potential as a function of 
the three auxiliary fields (a,ip T ,u T ), temperature T and chemical potential [i 



F cS =J£ ] + V q (m q ; fl, T) - N c log y^zZ , (2-74) 



+ \{Pt + 4/3 tsC 7 2 + 8/3 TT m'/)(^ _ ^) _ 4/3 TT m' q ip T a , (2-75) 



_ Ka + mp _ m' q 

m « ~ 1 + 4f3 TT <p T ' mq ~ JZ&l' [ ] 

Z± =1 + (p T + 2(3 TS a 2 + 8P TT m' )(<p T ± w r ) - A/3 TT m' q a , (2-77) 

6 CT =6' CT + 2/3 TT = 6 CT + 2/3 TT + 2/3 s a 2 + 6/3 ss a 4 + 2/3 rs (^ 2 - w 2 .) , (2-78) 



A =n - log . (2-79) 

Since tp T is a function of a and w r whose function form is independent from (T, /i) , 
we cannot regard it as an order parameter. Thus, we have two order parameters in 
the present treatment of NNLO SC-LQCD, the chiral condensate a and uj t . We 
can regard uj t is a vector potential field for quarks; the chemical potential shift is 
mainly determined by uj t , and uj t contributes repulsively to the effective potential in 
equilibrium. This two order parameter feature may be a natural consequence from 
the potential term from quarks, V q (m q ; fl, T). There are two independent derivatives, 
dV q /dm q and dV q /dfi, which appear in the equilibrium condition, then we have two 
degrees of freedom. 

§3. Chiral Phase Transition in NNLO SC-LQCD 

The effective potential derived in the previous section determines the equilibrium 
(vacuum) and the phase structure of QCD matter. In this section, we investigate 
the chiral phase transition in the chiral limit wiq = at iV c = 3. First, we discuss the 
effective potential surface and the stationary conditions. Then, critical temperature 
and chemical potential are investigated. Finally, we discuss the coupling dependence 
of the critical point. 

3.1. Effective potential surface 

The equilibrium is determined by the stationary condition of J- c r with re- 
spect to the auxiliary fields. By substituting ip T in Eq. (|2-74p with the solution of 
Eq. (|2-73p . we obtain the effective potential as a function of a and u T , F e ^{a,oj T ) = 
J-cs(a,u} T ,ip T = (p T (a,uj T )). The remaining stationary conditions are for a and u T 

dleS = <9Jeff = q / 3-1 > 

da du T 

and we have to solve this coupled equation self-consistently. This is equivalent to 
searching for the saddle point of J- e g EDJS)J33 where J- e g is convex downward and 
upward in a and uj t directions, d 2 J 7 e ff/da 2 > and d 2 T e s/doj 2 < 0. In Fig. [31 we 
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display T^s{(T, uj t ) at (T,/j,,(3) = (0.1, ^ c ,T=o.ii 5.0), as an example. The solid curve 
in Fig. [3] shows the stationary condition for u T , uj t = cof' at '(a), and the filled circles 
show the equilibrium points. The effective potential as a function of a for given 
(T, fx) is defined as F c s{a) = T e g(a,uj T = w^ tat '(<r)), whose minimum corresponds 
to the equilibrium. In the left (right) panel of Fig. [JJ we display J r e ff( cr ) on the 
T-(/i-)axis. The filled circles show the equilibrium points, which are obtained from 
the stationary condition for a, Eq. ()2-63p . 




Fig. 3. The effective potential as a function of a and uj t at (T, fj,, /3)=(0.1, fi c , 5.0) in the lattice 
unit. The stationary condition for ui T obtained in Eq. (|2-69|) is satisfied on the solid curve. 

In vacuum, the chiral condensate is non-zero, i.e. the chiral symmetry is sponta- 
neously broken and the equilibrium is in the NG phase. As temperature increases on 
the T-axis, the chiral condensate decreases smoothly and become zero at T = T CjA1= o, 
i.e. the chiral phase transition to the Wigner phase takes place. As chemical poten- 
tial increases on the /i-axis, the chiral condensate changes little and jump to zero at 
i 1 = Me,T=o- The chiral phase transitions are the second- and the first-order on the 
T- and the //-axes, respectively, at (3 = 5.0. These results are consistent with those 
in SCI j2D,[2!mi32j and NLC (34| SC-LQCD. 

3.2. Critical temperature and chemical potential 

The second-order phase transition boundary is obtained from the condition of 
C2 = where T c s(a) = ^2 n C n a n /nl. In order to obtain C2, we start from the 
effective potential as a function of a,<p T and uj t , JF e s(a,oj T ,ip T ) in Eq. (|2-74|) . From 
the stationary conditions, ip T and oj t are found to be even functions of a in the 
chiral limit (d<£>' /da\ a=0 = 0, = (p T ,u T ). Because of the stationary condition, the 
first derivative of J-" e ff by auxiliary fields are zero (dT e s/d^ = 0). By using these 
properties, C2 is expressed as, 

Id r 3^sV 

\ <P'=y T ,w T J 



8a 2 



(3-2) 
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Fig. 4. The effective potential in the lattice unit is illustrated as a function of the chiral condensate 
a on the T-axis (left) and the /x-axis (right) for {3 — 5.0. The filled circles show the equilibrium 
points. In the left (right) panel, the effective potential decreases for increasing T(fi). 



__ C (X) + c (v q ) _ 



a 



(^4 



a 



where, 



a 



(X) 



d 2 F 



cff 



da 2 



k+ 4/3 TS + 8/3 TT C£, (<p z T -u z T ) 



a 



(V 9 ) 



d 2 V q 
' da 2 



Pq 

2 



a 



(Z+) 



a 



z A 



(3-3) 

-VrCm' , (3-4) 

^1 N c (N c + l)(N c + 2) (C ml \ 2 

3T(2 cosh(N c jl/T) + N c + 1) ' \ Z x J ' 

(3-5) 

(3-6) 

In Eqs. (|3-3|) - (|3-7|) . the right hand sides have to be evaluated in the conditions a = 
and p q = —dVq/dfi. 

The second order critical temperature at fi = is obtained from the condition 
C2 = at fjL = as 



o2 <v 

6 CT 



3 TT C m >^ , 



C 2 

T n _ T (SCL) . 



(X) iV c G 



^4 



(3-8) 



■,(SCL) 



where T" c w ^ ; = d(N c + l)(N c + 2)/[6(iV c + 3)] denotes the critical temperature at 
ji = in the strong coupling limit. Note that at fj, = 0, ui T becomes zero (i.e. 
jl = 0). Therefore, the coefficients in T CjA1= o are the value at u T = 0, and then 



a 



decreases as f3 increases because of finite 



we find Z + = Z_ (C 2 
coupling effect. 

In Fig. [5l we compare the critical temperature in NLO and NNLO. In both NLO 
and NNLO, T Cj(Ji= q decreases as f3 increases, and T c ^ = q in NLO and NNLO exhibit 
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Fig. 5. The p dependences of the critical temperature at fi = (T CiM= o) in the lattice unit. We 
compare the NLO (the dashed line with circles) and NNLO (the solid line) results. The trian- 
gles represent the results of the critical temperature (T c ,^=o, open triangle) and the critical 
coupling (/3 C , filled triangles) obtained in Monte-Carlo simulations with one species of un- 
rooted staggered fermion: From the left, T C)M= o in the SCL with MDP simulations, 1 ^ /3 C at 
(N T ,m ) = (2,0.025)F (2,0.05),P (4,0.0)P (4,0.05),™™ and (8,0.0).™ 
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Fig. 6. The ft dependences of the critical chemical potential at T — (/i c ,T=o) in the lattice unit. 
In the left panel, we show NNLO results for the critical chemical potential (the solid line), the 
first-order critical chemical potential (the dashed line with crosses) and the second-order critical 
chemical potential (the dotted line with squares). In the right panel, we compare the NLO (the 
dashed line with circles) and NNLO (the solid line) results. 



almost the same values. We also show the MC results with one species of unrooted 
staggered fermion on the critical coupling (/3 C ) for given N T = 1/T = 2,4 and 8 at 
[i = (filled triangles) and the critical temperature in SCL: From the left, we show 
T c ,u=o m the SCL with monomer-dimer-polymer (MDP) simulations ,® j3 c at N T = 2 
(the quark mass m = 0.025),® N T = 2 (m = 0.05)® 8 3 x 4 lattice (with a chiral 
extrapolation),® 8 3 x 4 lattice (m = 0.05),® 6 3 x 4 lattice (m = 0.05),® and 
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N T = 8 (with a chiral extrapolation).^ The reduction of T c ^ = q is not enough to 
explain the MC results. In this formulation, we do not consider the deconfmement 
phase transition and the quarks are confined. T CjA1= o could decrease than NNLO 
results from the effects of the Polyakov loop which is an order parameter of the 
deconfmement phase transitionPSJ In NNLO, the phase transition at fi = stays 
to be the second order in the region j3 < 6. This point differs from the numerical 
simulation results, where the phase transition at \i = is shown to be the first order 
in the continuum regionPS The relation with the critical point position will be 
discussed in the next subsection. 




Fig. 7. In the left panel, we show the critical chemical potential at T = (/ii ls *' , the dashed line 
with circles). The result obtained by "approximate expression", (E ( q vac) +Sfi w )/2, is also shown 
for a comparison (the solid line). In the right panel, we compare at T — a mass modification 
Eq Vac ^ (the upper solid line) with a chemical potential modification Sfj, w (the lower dashed line) 
as a function of f3. In both panels, we use the lattice unit. 

In Fig. [6l the phase transition on ^t-axis is numerically found to be the first-order 
for P < 5.5. We find that ^t=o 1S no * l ar g er y affected by finite coupling effects. 
This is understood as follows. Taking the limit T — > in Eq. (|2-52j) . one can easily 
derive that the quark free energy reduces to 

_ f -N c E q {E q > p) 
Vq ~\ -N c fi {E q <p) ■ 

As explained after Eq. (|2-54p . E q represents the quark excitation energy. In this 
formulation, we do not consider the deconfmement phase transition and the quarks 
are confined. Then, the one and two quark excitation contribution disappears af- 
ter temporal link (Uq) integration. As a result, the factor N c accompanies E q in 
Eq. (|3-9|) . This means that the contribution from three co- moving quarks appear as 
a baryon in the effective potential at low T. 

At fi = (ji!£t=oi effective potentials at two local minima are the same, then the 
following relation holds 

~N C E™ + (-^iV c < G + jg } ' NG - N c log Z X NG 

- N 5u w + (- AT X - W 1 ^ (x) ' w w 
2 c p V 2 



N c n + -N c 5^ + ( -N c 5^ + - N c log ) , (3-10) 
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where the lhs and rhs show the effective potentiais in the Nambu-Goldstone (NG) 
and the Wigner phases, respectively. We numerically find that the lhs and rhs in the 
brackets have almost the same values. Then, the above relation is approximately 
represented as 

- l -N c E^ G ~ -Net + \n c 8^ , (3-11) 

These approximate expressions are explained based on the two-body interaction 
dominance, where the total potential energy amounts to be half of the single particle 
potential. This relation is exact when the effective potential is a quadratic function 
of a or oj t . In the NG phase, the quark number density is very small at T = 0, then 
the main contribution to the potential energy comes from a. In the Wigner phase, 
the chiral symmetry is restored, then the effective potential becomes a function of u T . 
Since there are contributions from other auxiliary fields and the quark free energy is 
not a linear function of a or u) T , the above relation is an approximate one. In the NG 
phase at T = 0, the quark number density is almost zero, and the quark excitation 
energy is almost the same as that in vacuum, E^ G (fi > 0) ~ E^ G (fi = 0) = E q yac \ 
In addition, we obtain uj t = p q ~ N c in the Wigner phase at T = 0. By using these 
facts, the Eq. (|3T1|) reduces to, 

^(lst,App.) = 1 J E (va C ) + ^Wj f (3 . 12) 

«5/i W = Sfi(a = 0,lo t = N c ) . (3-13) 

In the left panel of Fig. we compare the first order critical chemical potential with 
its approximate expression in Eq. (|3-12|) . For /3 > 3, fij^^ p '^ can roughly explain 

(1st) 

In the right panel of Fig. [7J we show the /3 dependence of E q and Sfj,. With 
increasing /3, E q becomes smaller but <5/U W goes to a larger value. The sum of them 
appearing in Eq. (|3-12p slightly increases and is not very sensitive to /3. Thus, the 
quark mass (or excitation energy) and chemical potential suppressions which come 
from the finite f3 effects cancel each other out, and we observe a small modification 
of the critical chemical potential. 

3.3. Critical point evolution 

In Fig. [8l we show the critical point (CP) evolution with (3 in NNLO (the solid 
line with circles) and NLO (the dashed line with circles). We find that the CP stays 
to be the tri-critical point (TCP) in NNLO. This point is different from NLO results, 
where TCP starts to deviate from the second order boundary at /3 ~ 4.5 and becomes 
critical end point (CEP) for larger ^jMUSS) While the temperature of CP decreases 
in both NNLO and NLO, we find different behavior in the chemical potential of CP, 
//CP- In NNLO (NLO), CP moves in the smaller (larger) \x direction with increasing 
(3 in the range > 3. The behavior in NNLO is favorable, since the phase transition 
with one species of unrooted staggered fermion at fi = is numerically shown to be 
the first order in the continuum regionPO 
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Fig. 8. The critical point (CP) flow resulting from increasing /? in T — ^1 plane in the lattice unit. 
The solid and the dotted lines represent trajectories of NNLO and NLO cases, respectively. We 
also show the first and second-order phase transition boundaries obtained in cases of SCL and 
NNLO at P = 6. 
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Fig. 9. The critical point flow with increasing j3 in several truncation schemes in the lattice unit. 
From the left panel, a part of NNLO effects, NNLO-TT2, NNLO-SS2, NNLO-TS12, NNLO-TT1 
or NNLO-SS1 is taken into account. We compare them with NLO (the dotted line with circles) 
and full NNLO (the solid line with circles) results. 



The above behavior of /icp is caused by the two temporal plaquette diagram, 
NNLO-TT2. In NNLO, the diagrams with connected two plaquette configurations 
shown in Fig.[2]contribute to the effective action. In Fig.[9l we show the effects of each 
NNLO diagram on the CP evolution. From the left, we consider the effective action 
terms only from the NNLO-TT2, NNLO-SS2, NNLO-TS1+NNLO-TS2, NNLO-TT1 
or NNLO-SS1 diagram. We find only the NNLO-TT2 shifts CP in the lower \i 
direction. The other NNLO diagrams do not modify the direction of CP evolution 
with (3. The temporal hopping of quarks, that is the thermal effect of quarks, seems 
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Fig. 10. The NNLO phase diagram at P = 6.0 in the lattice unit. The solid and thick dashed lines 
represent the first- and second-order phase transition boundaries, respectively. The filled circle 
represents the tri-critical point. On the (thick and thin) dashed line, C2 = is satisfied, where 
C2 is the curvature of the effective potential in Eq. (|3-2[) . Note that there is no phase transition 
on the thin dashed line between CP and X. By contrast, in the region where T is lower than 
that of X, two sequential phase transitions occur as fi becomes large. 



to be essential for the suppression of //cp a t larger f3. 

In Fig. [TUl we display the phase diagram at f3 = 6.0. The NG phase {a 7^ 
0) appears in low T and low fi region, and the Wigner phase (a = 0) appears 
at high T or high fi. We find that the second-order phase boundary obtained by 
C2 = exceeds the first-order one in low T region. Between those boundaries, 
we observe an interesting matter state where the effective chemical potential fi is 
adjusted to be around the quark excitation energy E q . In other words, the chiral 
symmetry is partially restored to satisfy the balance relation \x ~ E q , which leads to 
an intermediate value of the quark number density p q . In this partially chiral restored 
(PCR) matter, N c quarks are co-moving and they have a smaller constituent quark 
mass; the former shows the confined nature and the latter is due to the partial chiral 
restoration. When the PCR matter appears, two sequential phase transitions take 
place as fi becomes large with the balance relations jl ~ E q {^ rh q ). This phase 
transition pattern is consistent with the implication from the large iV c argument, 
i.e. the quarkyonic transition^ as discussed in our previous works.®'^ In order 
to clarify the relation of the PCR phase and the quarkyonic matter, it would be 
necessary to take account of the deconfinement transition, which will be reported 
elsewhere p3 

While it is still uncertain that the PCR matter realizes in the continuum limit, we 
emphasize that this third phase (besides NG and Wigner phases) has been commonly 
observed in the T-fi-fS space in the NLO and NNLO SC-LQCD via the self- 

consistent treatment of the vector potential u T . Compared to NLO resultspS'ED 1 
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PCR phase appears in the lower T region in the present work. NNLO-TT1 and 
NNLO-TT2 are found to be responsible for such a tendency. 

§4. Summary 

In this paper, we have derived an analytic expression of the effective potential 
at finite temperature and density including the next-to-next-to-leading (NNLO) ef- 
fects in the strong coupling expansion of the lattice QCD, and investigated NNLO 
effects on the chiral phase transition and its phase diagram. We adopt one species of 
unrooted staggered fermion corresponding to Nf = 4 in the continuum region. Ef- 
fective action terms have been systematically evaluated based on the strong coupling 
expansion. We have concentrated on the leading order of the 1/d expansion, then the 
NNLO effective action terms are generated from six types of two plaquette configu- 
rations. We have applied the extended Hubbard-Stratonovich transformation®'® 
in order to bosonize fermion interaction terms. We encounter those terms containing 
the interaction between the next-to-nearest neighbor sites (x, x + 20) arising from 
l/g 4 terms of the effective action. These effects can be evaluated by introducing 
the gluonic dressed fermion, which leads to modifications of temporal quark hop- 
ping and mass terms. We have obtained the effective potential as a function of 
temperature (T), chemical potential (fi) and the two order parameters: the chiral 
condensate (a) and the vector potential (oj t ). The equilibrium is determined from 
the stationary condition of the effective potential with respect to the auxiliary fields. 
NLO and NNLO effects result in modification of the wave function renormalization 
factor, quark mass and chemical potential. The effective chemical potential reflects 
the repulsive contribution of the vector potential. 

We have found that the critical temperature at // = (T Ci/J= o) is largely sup- 
pressed with increasing /?. In comparison, the chemical potential at T = (/i Cj T=o) 
is a slightly increasing function of /3. These behaviors of the critical temperature and 
chemical potential in NNLO are consistent with those in NLO in the region (3 < 6.0. 
The behaviors h c ,t=o an d T c ^ = q are understood from the quark mass reduction, and 
its cancellation with quark chemical potential reduction, respectively. 

We have found that the critical point (CP) moves to lower T direction with 
increasing (3 in both NLO and NNLO. In NNLO, the chemical potential of the 
critical point decreases with increasing f3. This shift is opposite to the NLO case 
with (3 > 3.0.®'® Thus, the critical point (CP) flow with increasing f3 is sensitive 
to l/g 4 effects, and the first-order transition tends to dominate the phase boundary 
as P is increased in NNLO. The NNLO CP flow would be favorable because the 
phase transition at fi = is shown to be the first order for the unrooted staggered 
fermion in the continuum region.® We have found that the next-to-nearest neighbor 
interaction makes the CP flow go to a lower /i direction. 

In low T and high /i region, we have found the partially chiral restored (PCR) 
phase where effective chemical potential is always adjusted to be around the quark 
excitation energy. The PCR matter is obtained also in NLO as long as the two order 
parameter (a and u T ) are introduced. The appearance of PCR matter would be a 
common consequence in NLO and NNLO SC-LQCD with finite coupling effects. 
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There are several points to be studied in future. Firstly, the Polyakov loop can 
be constructed by using two plaquettes in the NNLO SC-LQCD with iVV = 2. This 
effect has not been considered in the present analyses, and quarks are confined. It is 
interesting to investigate the deconfinement transition in addition to the chiral phase 
transition by considering the Polyakov loop effects. Secondly, we should study the 
higher order of the \ jd expansion, which contains the spatial baryon hoppings and 
could play an essential role at finite density on the lattice. 
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